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ABSTRACT 


This thesis addresses the optimization of all types of space electrical propulsion 
thrusters. From the Langmuir-Irving payload mass fraction formulation, a “dual- | 
optimum” solution is defined, yielding a minimum overall mass for a specified payload 
consistent with minimum transfer time. This solution fixes the ideal payload mass ratio 
(mp / Mo) at a value of 0.45, establishing the ratios of effective exhaust velocity (v / v<) 
and incremental change of vehicle velocity (Au / v,) to characteristic velocity at 0.820 
and 0.327 respectively. The characteristic velocity (v<) includes thrust time as well as 
engine efficiency (n,) and specific power (a). A range of mass ratios from 0.35 to 0.55 is 
used in order to allow the system designer some flexibility while remaining close to 
optimal. Nine examples are presented which demonstrate that mission profiles can be 
optimized by profile-to-thruster matching. A comprehensive list of currently available 
electric propulsion engines is provided. This list details important parameters such as the 
specific power, which “sizes” an engine in terms of power provided to the thruster at the 
cost of additional mass. Allowance is also made for a fuel tank mass penalty, and 


examples show that this can also noticeably influence the optimum design. 
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I. INTRODUCTION 


A. SPACE ELECTRIC PROPULSION COMES OF AGE 

The application of electric propulsion for space is not a recent revolution in the 
field. The idea was proposed as early as the 1950’s. One significant article appeared as a 
chapter in the text, Space Technology, entitled “Low Thrust Flight: Constant Exhaust 
Velocity in Field-free Space,” written by D. B. Langmuir (Langmuir, 1959, pg. 9-01). 
Aside from the evolution of technology, little has changed in the fundamentals of electric 
propulsion since that time. 

In the last 33 years, more than 300 electric propulsion thrusters have flown on 
over 100 spacecraft, and in the next decade, a significant increase in electrical propulsion 
employment is expected (Filliben, 1997, pg. 5). Nonetheless, electrical propulsion is 
only now seeing rapid and wider-scale introduction into spacecraft propulsion programs. 

The first electrical propulsion engines flew in the mid-1960’s. The former Soviet 
Union used a pulsed plasma thruster (PPT) to provide attitude control to their Zond-2/1 
satellite, which launched in 1964 (Martinez-Sanchez and Pollard, 1998, pg. 692). Then, 
in 1965, the United States followed up their effort using a resistojet engine for orbit 
adjustment of the Vela/2 satellite (Martinez-Sanchez and Pollard, 1998, pg. 689). 

If both of those missions proved successful, why has additional implementation of 
electrical propulsion for space applications — so seemingly slow, spanning almost 40 
years? For this, there are both technical and ace reasons. 

First, engineers have only recently — able to package sufficient electrical 


power and power conditioning equipment onboard spacecraft in order to provide for the 





demands of the electrical propulsion system. Unlike chemical propulsion, electrical 
requires energy accumulation and energy processing functions creating a power supply 
penalty at the expense of payload mass (Humble, Henry, and Larson, 1995, pg. 510). 
Likewise, many mission planners have been reluctant to replace proven chemical rockets 
with comparatively untested electrical systems. As a result, electrical systems have only 
gradually become accepted and only for less demanding tasks -- such as north-south 
station keeping (NSSK) and orbit insertion. 

Nonetheless, electrical propulsion has several advantages over heentes rockets 
for near-zero gravity applications. Moreover, for deep-space fica. which require 
huge incremental changes in vehicle velocity (the factor Au), electrical propulsion is the 
only realistic option. In fact, it is likely that the future of deep space exploration will 
depend largely upon siectiic propulsion systems as chemical systems will simply be 


unable to attain the high Au required for these missions. Table 1 (Typical Mission 


Velocity Requirements) presents a summary of various mission Au requirements. 


Earth to LEO’ 7,600 m/s . 
LEOtoGEO” 4,200 
LEO to Earth escape 3,200 
LEO to lunar orbit (7 days) 3,900 
LEO to Mars orbit (0.7 yr) 5,700 
LEO to Mars orbit (40 days) 85,000 


LEO to Neptune orbit (29.9 yr) 13,400 
LEO to Neptune orbit (50 yr) 70,000 
LEO to solar escape 8,700 
LEO to 1000 AU* (50 yr) 142,000 
LEO to Alpha-Centauri (50 yr) 30x 10° 


* Low Earth Orbit = 270 km 
° Geosynchronous Earth Orbit = 42,227 km radius 
° Astronomical Units (1 AU = 149,558,000 km) 


Table 1: Typical Mission Velocity Requirements. From Hill, Peterson, page 508. 











Note however, that these mission velocity requirements are calculated for chemical 
propulsion and not electric propulsion systems. Regardless, Table ] illustrates the typical 
magnitude of Au required for various missions as the numbers do not change 
significantly. The values will be slightly higher for electrical systems. 

On the down side, electrical propulsion systems cannot sveivonie high-gravity 
fields due to their low acceleration or low thrust. Therefore, the initial boost to attain 
escape velocity from Earth will need to be accomplished via chemical rockets, and it is 
evident that the types of missions slated for the future will be accomplished via a 
combination of both electrical and chemical propulsion. There are specific performance 
limitations to each. (Hill and Peterson, 1992, pe. 491) 

Chemical rockets can be described as “energy limited” in that the fundamental 
chemical behavior of the propellant defines the quantity of energy (per unit mass of 
propellant) that can be released during combustion. High propellant energy for 
spacecraft is possible if a separate energy source is utilized -- such as nuclear or solar 
energy. In that manner, the rate of conversion of nuclear or solar energy to electrical and 
then to propellant kinetic energy is limited by the mass of conversion equipment required. 
That mass can be a comparatively large portion of the total vehicle mass. As a result, the 
electrical rocket is generally said to be “power limited.” (Hill and Peterson, 1992, 
pg. 491) 

B. PHYSICS BEHIND ELECTRIC PROPULSION 

Space electrical propulsion systems consist of three basic subsystems. Those 

subsystems include the power processing unit, sronellant management assembly, and the 


thruster assembly consisting of one or more thrusters. 


The power processing unit (PPU) contains electrical interfaces that link the 
spacecraft power source and the propulsion system. The PPU sends power and command 
signals to the thruster(s) and propellant flow control signals to the propellant subsystem. 
Additionally, the PPU is responsible for converting various power sources (solar, 
batteries, fuel site nuclear, etc.) to the proper voltage, frequency, pulse rate, and current 
suitable for the system components. 

The propellant management assembly (PMA) provides the means of storing, 
metering and delivering the propellant. The thrusters convert electrical energy to kinetic 
energy of a propellant in the form of exhaust or thrust. 

Regardless of type, the common feature of electric propulsion systems is the 
addition of energy to a working fluid from an electrical source. Operating in a steady or 
pulsed mode, the exhaust is accelerated via one or more of the three fundamental types of 
thrusters. Those types of thrusters are (1) electrothermal, (2) electrostatic, and (3) 
electromagnetic. (Martinez-Sanchez and Pollard, 1998, pg. 688) 

In reading the following sections, reference Figure 1, “Operating Principles of 
Various Electrical Propulsion Thrusters,” found at the end of this chapter on page 14. 
The figure will prove exceptionally useful in understanding the operation of several of 
the following thrusters as they are | 

1. Electrothermal! Thrusters 

In the electrothermal thruster, the propellant is electrically heated and 
thermodynamically expanded. The heated gas then achieves supersonic velocity as it 


passes through a nozzle. Overall, these thrusters have one major drawback in that, as a 


class, they produce limited high exhaust velocities. Additionally, their material 











characteristics limit their performance to values that are generally in line with those 
exhibited by chemical rockets. At present, there are two basic types of electrothermal 
see -- resistojets and arcjets. 

The resistojet is the more elementary of the two in its class. In essence, it 
operates by using an electric current to directly heat metal components (coiled wire, 
tubing, or fins) which then transfer heat to the gaseous propellant through radiation 
and/or convection. As a result, the resistojet’s performance is limited by the capacity of 
its structural components to withstand high temperatures, and resistojets achieve only a 
modest specific impulse of 300-310 seconds. This is due to the relatively high molecular 
mass of the propellant gas in addition to material limitations of the heating wall, which 
can sustain temperatures up to approximately 2000 K. (Martinez-Sanchez and Pollard, 
1998, pg. 689) 

Even with its comparatively lower value of specific impulse, the resistojet’s 
superior efficiency contributes to far higher values of thrust/power than other varieties of 
electrical propulsion. Additionally, these engines possess the lowest overall system “dry 
mass” since they do not require a power processor and their plumes are uncharged. 
(Sackheim and Byers, 1998, pg. 672) 

Resistojets have seen successful employment on such craft as Intelsat V, Satcom 
1-R, GOMS, Meteor 3-1, Gstar-3, and Iridium S/C in addition to older satellites and test 
flights (Martinez-Sanchez and Pollard, 1998, pg. 689). They are most attractive to the | 
mission planner for low-to-modest energy applications. In particular, they are most 
suitable where power limits, thrusting times, and plume impacts are mission drivers. 


(Sackheim and Byers, 1998, pg. 672) 





The other form of the electrothermal thruster, the arcjet, manages to overcome the 
thermal problems encountered by the resistojet, and it enables a significant increase 
(approximately double) in specific impulse over the resistojet. The arcjet accomplishes 
this by forming an attached arc within the nozzle that directly heats the propellant stream 
to temperatures much greater than those encountered at the thruster body. Essentially, 
the power is deposited internally in the form of the electric arc, which typically operates 
between a concentric upstream rod cathode and a downstream anode that serves as a 
supersonic nozzle. (Martinez-Sanchez and Pollard, 1998, pg. 689) 

The arc core reaches temperatures of 10-20,000 K while the buffer layer near the 
wall is maintained at less than 2000 K. As a consequence of the temperature gradient, 
there is practically no propellant flow through the are core. This results in high specific 
impulse at the cost of reduced propulsive efficiency. In addition, the flow structure at the 
throat is extremely non-uniform, and the thruster is primarily limited by arc instabilities 
and erosion of the nozzle. (Martinez-Sanchez and Pollard, 1998, pg. 689) Generally, 
arcjets exhibit significant decreases (approximately six times) in thrust/power relative to 
the resistojet. This is due to an increased specific impulse coupled with relatively low 
values of efficiency ranging from 0.3 to 0.4. (Sackheim and Byers, 1998, pg. 672) 

One significant disadvantage to the arcjet is that the required PPU is substantially 
more complex than that encountered in the fesisivie: as the complex plasma arc must be 
controlled. In fact, the PPU may be several times heavier than the thruster itself: Still, 
however, the arcjets exhibit relatively low dry masses when compared with some of the 


other electrical propulsion systems in use or development. (Martinez-Sanchez and 


Pollard, 1998, pg. 689) 








With its intermediate specific impulse of 600 - 650 seconds, the arcjet is a very 
viable option for short burn duration missions. The lower specific impulse implies higher 
thrust for a given power, and arcjets have been used for geostationary applications, 
including the Telstar 4 and GE-1 satellites. (Martinez-Sanchez and Pollard, 1998, pg. 
689) The engines are relatively simple to integrate, and they are the least costly and 
complex of any plasma propulsion device (Sackheim and Byers, 1998, pg. 672). 

2. Electrostatic Thrusters 

Electrostatic thrusters rely on electrostatic fields in order to exert forces on the 
propellant which is throttled into the field in the form of charged particles. Such thrusters 
can only operate in a vacuum, and all the created particles must be of the same sign -- 1.e. 
all positive or all negative. Likewise, those particles must be neutralized after passing 
through the engine. Otherwise, the thruster would eventually render itself ineffective as a 
net charge would build up on the spacecraft setting up the return of oppositely charged 
~ ions to the craft. As a result, any thrust would be cancelled. In addition, sensitive system 
components could be'damaged by the returning ions. (Humble, Henry, and Larson, 1995, 
pg. 532) 

Naturally, electrons wail seem to be the most likely candidate particles for these 
engines. They are easy to produce and to accelerate; however, their small mass renders 
them impractical. They are, in fact, so extremely light that the momentum imparted to 
them is negligible even with their high velocity. 

Consequently, electrostatic thrusters operate using heavy atoms charged as 
positive ions or as colloids (liquid droplets). Neutralization of the positive ions is easily 


accomplished with electrons, and the ions can be 240,000 times as heavy as the electron. 








The charged colloids are generally at least 10,000 times as massive as the ions. Thus, 
substantial momentum can be generated. 

The implementation of heavy atoms has other advantages. The thrust per unit 
area increases as the square of the particle mass to the charge ratio in addition to the 
highly desirable characteristic of high voltage and low current as opposed to low voltage 
and high current (Sutton, 1992, pg. 580). | 

Electrostatic thrusters are categorized by the source of charged particles, and all 
three types exhibit the same three basic stages: ion production, acceleration, and © 
neutralization. The five types of thrusters that have been studied are as follows: 

1. Electron Bombardment Ion Thruster. Positive ions are produced via the 
bombardment of vaporized or gaseous propellant (such as xenon or mercury) 
with electrons emitted from a heated cathode. 

2. Ion Contact Thruster. Positive ions are created by passing a propellant vapor 
(cesium) through a hot contact ionizer (usually tungsten). Although it was 
intensely navestiosied 15 to 30 years ago, the ion contact thruster was 
abandoned as impractical. 

3. Field Emission (Colloid) Thruster. This engine generates positive or 
negative particles by passing tiny liquid droplets of propellant through an 
intense electric field (corona discharge). (Sutton, 1992, pg. 580) 

4. Radio Frequency Ion Thruster. This thruster produces a high specific impulse 
at the price of saplenity and mass. Besentally, it consists of a discharge 


chamber that is surrounded by an inductive coil connected to a radio 


frequency generator. Neutral propellant atoms are ionized by the 











bombardment of electrons which are accelerated by the induced high- 
frequency eddy fields in the discharge chamber. 

5. Electron Cyclotron Resonance / Microwave Discharge Ion Thrusters. Here a 
circularly polarized microwave beam is projected from a conventional 
waveguide via a suitable dielectric window in order to ionize the propellant in 
the discharge chamber. (Filliben, 1997, pg. 5) 

At the present time, 1on engines are the “thruster of choice” for deep missions 
such as interplanetary transfers. These missions require a high Au and must tolerate long 
thrusting times. (Martinez-Sanchez and Pollard, 1998, pg. 692) A typical xenon ion 
thruster can generate an average mission specific impulse of 2,800 - 3,500 seconds with a 
lifetime of greater than 8000 hours (333 days). . 

3. Electromagnetic Thrusters 

There are several variations to the electromagnetic type of thruster. All make use 
of a propellant gas that is heated to a plasma state. That plasma is then part of a current- 
carrying electric circuit which interacts with combined electric and magnetic fields 
thereby generating thrust. 

The four types of electromagnetic thrusters which we will discuss include the 
Hall-effect thruster, stationary plasma thruster (SPT), the magnetoplasmadynamic (MPD) 
thruster, and the pulsed plasma thruster (PPT). 

The Hall-effect thruster is named for its closed circular electron drift that is 
exhibited between the cathode and anode of the thruster. This effect is analogous to the 


Hall Effect. The SPT operates on similar principles. 


In the Hall-effect thruster, an axial electric field and radial magnetic field are 
established in the discharge chamber. Meanwhile, neutral propellant atoms are fed into 
the discharge chamber and ionized by the bombardment of electrons. The siaieiias 
electric and magnetic fields then force the electrons, emitted from an external hollow 
cathode, to become trapped in azimuthal drift motion. At the same time, the radial 
magnetic field is not sufficiently strong to affect the propellant ion trajectory, and they 
are consequently accelerated scaly by the electric field, producing thrust. (Filliben, 
1997, pg. 5) 

- The first use of a Hall thruster electric propulsion system in space took place in 
October of 1998. The U. S. Navy’s Hall thruster, installed on an experimental spacecraft 
owned by the National Reconnaissance Office, began operating as part of the electric 
propulsion demonstration module (EPDM). The engine provides for orbit raising and 
long-term station keeping. (Space Tracks, 1999, pg. 13) 

The particular engine onboard the EPDM is a TAL (thruster with anode layer) 
type Hall thruster, using xenon gas. It is a gridless device producing greater than twice 
the thrust of an ion engine of the same operating power. (Space Tracks, 1999, pg. 13) 

Although still in ‘eecioaiient the magnetoplasmadynamic thruster (MPD) is 
generally regarded as the leading potential candidate for future deep space missions such 
as , heavy-lift Mars transfer, utilizing a nuclear powerplant. (Martinez-Sanchez and 
Pollard, 1998, pg. 693) 

The engine operates by generating a current along a conducting bar. The current 


consequently yields a self-induced azimuthal magnetic field that interacts with the current 


of an arc that travels from the point of the bar to a conducting wall. The resulting force 











has two components. One component is a radially inward force that constricts the flow, | 
and the other force acts along the axis, producing directed thrust. As expected, one of the 
great problems with this thruster is the substantial seen that occurs at the point of 
contact between the arc current and electrodes. (Santarius, 1997) Additionally, the MPD 
thruster is more difficult to optimize due to the fact that several physical effects are 
involved in its operation. 

The Pulsed Plasma Thruster (PPT) is “different from all other concepts” in two 
fundamental ways. This device operates in short pulses of approximately 10 
microseconds duration, and in its most developed form uses a portion of solid propellant 
feedstock (typically Teflon®). As a result, no propellant tanks (referred to as “tankage”) 
are required in this device. (Martinez-Sanchez and Pollard, 1998, pg. 692) 

In the PPT, a portion of the feedstock is ablated and ionized by an electrical arc 
discharge sheet initiated between two electrodes by a discharging capacitor. The 
resulting propellant plasma is accelerated by the interaction of the arc and the self- 
induced magnetic field of the current loop. (Filliben, 1997, pg. 5) 

The PPT design is comparatively simple and reliable. However, one of the 
inherent problems with this thruster is that the force falls off as the current loops get 
large. In addition, PPT’s generally suffer from low efficiencies on the order of 8 to 13 
percent. (Santarius, 1997) 

PPT’s offer two additional distinct advantages in that they integrate a non-toxic 
propellant feed system with a thruster in a single compact unit as well as offering a — 
variable pulse repetition rate. As a result, they possess a flexibility of operation over a 


wide range of mean power or thrust, and therefore, they have found employment for 


1] 


precision orbital or attitude-control lai since the late 1960’s with the LES-6 satellite 
and the U. S. Navy’s NOVA constellation. (Martinez-Sanchez and Pollard, 1998, 
~ pg. 692) 

The optimal implementation of these three fundamental types of electric 
propulsion engines (electrothermal, electrostatic, and electromagnetic) will be the focus 
of our analysis. As illustrated in this chapter, ihe are several variants of each of the 
engines. Some of those engines are currently in use OF in various stages of development. 
Please reference Table 2 on page 13 for a summary of typical electrical propulsion 


performance parameters. 
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Figure 1. Operating Principles of Various Electrical Propulsion Thrusters. 
From Martinez-Sanchez and Pollard, page 691. 


(a) resistojets 
(b) arcjets 

(c) Hall thrusters 

(d) ion engines 

(e) pulsed plasma thrusters 

(f) field-effect electrostatic propulsion thrusters 


(g) self-field magnetoplasmadynamic thrusters 
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iI. THEORETICAL BACKGROUND 


A. FIRST PRINCIPLES 

Electrical propulsion systems have four general yet unique characteristics. These 
characteristics are exploited in analyzing and understanding the electric propulsion 
system, and they are delineated in the following paragraph. 

Electric propulsion engines (1) incorporate a comparatively heavy mass of 
required power generation equipment, and they (2) produce low thrust. As a result of that 
low thrust, (3) very low accelerations (on the order of 10% to 10° g,) are attained (4) 
requiring long operating times. In the following development of the equations governing 
electrical engines, the electrical rocket’s performance will be chararienve in terms of 
power and mass as first proposed by Langmuir and Irving. (Sutton, 1992, pg. 596) 

First, the mass ante of electrical propulsion will be analyzed. 

The major components integral to an electrically powered space vehicle are 
broken down as follows: 

Mo = Mp + Mp) + -Mpp. (1) 
In words, this equation states that the initial total mass of the vehicle (mg) is equivalent to 
thie sum of the masses of the propellant (m,), payload (m,)), and powerplant (mpp). Refer 
to the list of symbols and abbreviations (on page xiii) for the comprehensive table of the 
following abbreviations and their respective ili 

The determination of mpp for an assigned mission is the most amorphous. This 
term incorporates a wide variety of components that make up the empty (or un-fueled) 


propulsion system. Those components include the thruster, the propulsion storage and 





feed system, the energy source with its conversion system and auxiliaries as well as the 
associated structure (Sutton, 1992, pg. 597). 

Now, let us look at the power aspects. 

Regardless of the energy source (battery, nuclear, solar, fuel cells, etc.) to the 
input power supply, the input energy must always be greater than the required electrical 
power output, and power conversion efficiencies typically average approximately 70 
percent. This aspect is canola for in the engine parameter known as specific power, 
a, or the “state-of-the-art.” It is defined as the ratio of electrical power output to the mass 
of the power plant, and it serves as i. indication of the system’s ability to effectively and 
efficiently generate or convert to electrical power. The equation for a is defined as 
follows: 

0. = Pe/ mpp, | | — Q) 
where P, is the electrical power that is actually supplied to the propulsion system. | 
Typical high “alpha” values that are currently attainable range from approximately 100 to 
200 W/kg. It is a singular goal in electrical propulsion to strive to push this number to 


the highest attainable values. (Sutton, 1992, pg. 597) 
In an electrical propulsion thruster, the electrical power input is converted to the 


kinetic energy of the propellant exhaust. Thus, P, can also be defined by accounting for 


losses due to thruster efficiency, n:, where 


m= Pye / Pe = (1/2) m v’/P, (Sutton, 1992, pg. 571), (3) 


and, consequently, we may redefine P, as the following: 


Pe =O Mpp = (1/2) m V"/ m= My V'/ (2 ty NH). (4) 











In this equation, m is the propellant mass flow rate, v is the effective exhaust velocity, 
and tp is the duration of operation in which the propellant mass, mp, is ejected from the 
thruster. It is important to note that this equation assumes a uniform rate of propellant 
expulsion. (Sutton, 1992, pg. 597) | 

Using the previous definitions of Mo, O&, and Pe, we can form the equation for the 
payload mass fraction, which is written, 

Mp} / Mo = (1 - (V*M(2 & ty H(A” - 1))/ 8”. (5) 
This equation assumes sean and drag-free flight. Additionally, by defining the 
characteristic velocity, v<, as 

V.=20 ty Nt (6) 
we can re-write the payload mass fraction as follows: 

mp / Mo = (1 = (v /ve?(e8™ - 1)y/ eb. re) 
In this form, the payload mass fraction is more manageable. (Sutton, 1992, pg. 597) 

A plot of this equation is illustrated in Figure 2. The equation is graphed over the 
range of v/v, for various values of ” /v,. It is evident that a maximum exists such that 
ee ae ee ratio, Mp! / Mo, for a given velocity increment, Au / vg, is maximum at a 
particular value of exhaust velocity, v / v,, or specific impulse, I; / vo. Specific impulse 
and effective exhaust velocity are related via the equation, 

I, = v/ > where g5= 9.81 m/ 57, (Sutton, 1992, pg. 597) (8) 

The maximum described in these curves exists for two reasons. First, the inert 
mass of the power plant, mpp, increases almost linearly with the specific impulse. 
However, the mass of consumed propellant decreases with specific impulse. The vehicle 


velocity increment, Au, increases with higher I,; however, Au is also decreased by the 





lower mass ratio due to the increased inert mass. Consequently, there is.a unique 


optimum value of I, (or v) for every electric propulsion mission as defined by the 


corresponding required value of Au. (Sutton, 1992, pg. 597) 
Again, refer to Table 1 for a brief summary of various mission Au values 
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Figure 2. Payload Mass Fraction Plot: mp) / mg vs. Vv / Ve. 


The payload mass ratio can also be re-written to a somewhat more useful form by 
(9) 


solving for Au/v,. In this case, the equation takes on the following form: 


Au/v_=(v/ ve) In [(1 + (v/ ve)’)/((mp) / mo) + (V/ Ve)")], 
and it is plotted in Figure 3. In this case, the equation is graphed over the same range of 
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v/v, for various values of mp / mo. 


. Mp /m, = 0.35 





Figure 3. Payload Mass Fraction Plot: A u/v_. vs. v/ Ve. 


This graph shows the optimum curves for three different mass payload ratios. 
From this plot and the defining equations, we note that the optimal values for Au / v, and 


v/v, are as follows: 


m,; / Ms, Au/v.{opt v/v,.| opt 
0.35 0.404 - 0.782 
0.45 0.327 0.820 
0.55 0.257 0.864 


Table 3. Optimum Values: Au / v, and v/v, for Various Payload Mass Ratios. 
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These values will be an essential aspect of our engine selection (or engine-mission 
matching) algorithm, which is developed in Chapter IV. 
B. ANALYSIS — DUAL OPTIMUM METHODOLOGY 

Incorporating the previous development via “first principles,” it is now 
advantageous to establish a “joint” or “dual” optimum, which will optimizes the payload 


mass ratio and yields the minimum transfer time. It can be calculated by multiplying the 


payload mass ratio and the previously defined equation for Au / Ve. 
The plot of the resulting equation (for various payload ratios) is illustrated in 


Figure 4, and it shows that the dual-optimum payload mass ratio peaks at a value of 


My! / My = 0.45. 
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Figure 4. Dual Optimum. 
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This plot indicates that the optimum payload ratio is 0.45. However, due to the 
fact that it may be difficult for the designer to achieve this mass ratio, an approximate 
design range of payload mass ratios of 0.35 to 0.55 is investigated. As illustrated in 


Figure 4, these values do not lie too far below the optimum for 0.45. 


It should be restated that for constant propellant mass flow (m ) constant 


thrust (F) and negligibly short start and stop transients: 


I,=F/m g. and «.F=m gol,. (10) 

Hence, for a higher specific impulse, I,, less fuel is required. However, as the specific 
impulse is increased, we may drive through the optimum due to the equation for Au/ v,, 
_ which incorporates the applicable “efficiencies” of the engines in the forms of «, tp, and 
nt As aresult, an engine with an arbitrarily high I, may not be optimum for a particular 
application or mission. For example, the Teflon PPT has relatively high values for I, at 
the expense of comparatively poor values for a and 1. 

Figure 5 illustrates a three-dimensional optimum surface. The plot is illustrated 


across the axis values of... 


X! V/Ve_ range: 0.5 to 1.3 
y: Mp1 / Mo range: 0.3 to 0.6 
Z: (Mpi / mMo)(Au / Vc) range: 0.115 to maximum at 0.148 


The plot clearly shows the establishment of a optimum or maximum point at the payload 


mass ratio of precisely mpi / Mp = 0.45. 
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Figure 5. Three-dimensional Dual Optimum Surface. 


C. TANKAGE PENALTY 

The parametric curves that have been presented up to this point are inherently 
ideal in nature. They generally assume that all spacecraft mass can fall neatly into the 
three categories of mass of propellant, payload, or power-plant. However, as previously 
noted, it can be difficult to assess where exactly certain masses (such as shared power 


conditioning equipment) should be charged, and as a result, the calculation of the specific 


power, a, for the respective electric propulsion engine can have high variability. In the 
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same manner, one can assess somewhat of a “tankage penalty” in proportion to the 
amount of propellant required as a result of the following reasoning. 

The typical propellant tanks are composed of overwrapped titanium or aluminum. 
Their loading varies depending upon several factors including satellite mass, expected 
life on orbit, mission (orbit raising, NSSK, etc.) and cant angle of thruster. As a result, 
the actual size and mass of tanks will vary regardless of engine type or manufacturer. 

An accurate figure of merit for the required tankage revolves around the following 
assumptions. First, that a propellant tank that holds zero propellant mass is assessed at 
one kilogram. Second, that the tank has a mass fraction of approximately ten percent as a 
function of loaded propellant. It 1s important to note that all of the listed electrical 
propulsion engines, except the pulsed MPD type, require tankage. (Clauss, 1999) 

For example, a tank that holds 40 kg of xenon would have a dry mass of... 
1.0kg+ 0.1 (40kg)=5 kg. | 

If this 10 % tankage “penalty” and one kilogram of dead dry mass are 
incorporated into the previous analysis, the equation for the initial mass and components 
is altered from Mp =Mp + Mp! + Mpp to Mo= 1.1 Mp +My, + Mp). However, the 
equation for final mass remains unaltered as follows: mg=mp- mp. The 0.1 mp is 
residual mass that remains within the system. 

Using these definitions, the equation for the rewritten payload mass ratio, solved 
for Au / v,, takes the following form: 

Au / Vo = (v/ Vc) In [(1.1 + (v / vc)?/((mpt / mo) + (v / Ve)” + 0.19]. (11) 


This equation is plotted as Au / v, vs. v/v, in Figure 6. 
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Figure 6. Payload Mass Fraction Plot with 10 % Tankage Penalty: 


Au/ V¢_ vs. V/ Ve. 


Note that the curve is essentially similar to the plot illustrated in Figure 3; 
however, it is adjusted for the ben percent tankage “penalty.” As a result, the optimum 
curves have shifted down and right from their original values. This is illustrated in 
Figure 7 which shows the optimum curves with and without the tankage “penalty” for the 
payload ae ratio (Mp; / Mp) of 0.45. 


For a ten percent tankage “penalty,” the following optimal values for Au / v, and 


v/v,-result: 
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Mp\/My Au/v,jopt v/vc-| opt 


0.35 0.375 0.850 
0.45 0.306 0.895 
0.55 0.240 0.945 


Table 4. Optimal Values with 10% Tankage Penalty: Au/v, and v/ ve. 


My /m,=0.45 


| / : Mpi /m, = 0.45 with 10% tankage penalty . 
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Figure 7. Optimum Comparison with and without 10 % Tankage Penalty: 


Au / Vv, VS. V/ Ve at Mpi / Mo = 0.45. 


For all of the plotted “tankage” penalty curves, the values plotted are fora l0% © 
penalty. Typical tankage fractions range from 0.05 for the hydrazine (N2H4) resistojet 
and arcjet to 0.15 for the Hz arcjet. Xenon Hall thrusters and ion thrusters have tankage 


fractions of about 0.12. (Martinez-Sanchez and Pollard, 1998, pg. 690) 
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Il. CURRENT TECHNOLOGY -- ELECTRIC ENGINES 


A. NASA DEVELOPMENT AND OBJECTIVES 

Today, there are many organizations pursuing the research and development of 
space electric propulsion. Various private companies (including Primex, Hughes, ARC, 
and Daimler Chrysler) presently offer electric propulsion engines for sale and use in 
various applications. Additionally, several research facilities -- such as NRL (Naval 
Research Labs) and NASA’s JPL (Jet Propulsion Labs) -- have also developed various 
engines or components of engines. However, the one department of NASA that is 
primarily responsible for electric propulsion research is the On-board Propulsion Branch. 

NASA’s On-board Propulsion Branch is based at the John H. Glenn Research 
Center located on Lewis Field in Cleveland, Ohio. The branch is responsible for 
developing advanced on-board spacecraft propulsion technologies -- both chemical and 
electrical -- for future NASA missions. The electrical propulsion technologies currently 
under investigation at Lewis Field include resistojets, arcjets, gridded ion xenon thrusters, 
Hall xenon ion thrusters, and pulsed rn thrusters (PPT’s). (Oleson, 1999) 

In October of 1998, the on-board propulsion branch celebrated two major 
milestones in the realm of electrical eeopulcion: Deep Space One (DS1) was successfully 
launched on October 24, 1998 from Florida’s Cape Canaveral Air Station. This craft is | 
hall to validate new space technologies, including its ion propulsion primary engine, 
which was developed in the NSTAR (NASA Solar Electric Propulsion Technology 
Application Readiness) program. NASA Glenn developed the engine and power 


processors for DS1, and this flight was the first-ever application of such a technology to a 


2 


deep space mission. DS1 was to rendezvous with the near-Earth asteroid, 1992 KD, in 
July 1999, and it passed within 15 miles of the asteroid Braille on July 28 at a speed of 
greater than 35,000 mph (15.5 km/s). (Oleson, 1999) 

DS °s NSTAR ii engine utilizes 81 kilograms of xenon propellant, which it 
ionizes and accelerates to 68,000 mph. The engine, which is 12 inches in diameter, can 
operate at half throttle for greater than 20 months, generating a continuous thrust of 
0.02 lbf. Asa result, the engine is ten-times more fuel efficient than a comparable 
chemical system, and it is the prototype for smaller / lower-cost space vehicles. (Oleson, 
1999) 

On October 3, 1998, a STEX spacecraft was launched, carrying the EPDM 
(Electric Propulsion Demonstration Module). This coli was developed under joint 
effort on the part of NASA, Naval Research Labs (NRL) and private industry. The Glenn 
Research Center led the peseran, which was intended - demonstrate the first use of a 
high specific impulse, low-power Hall thruster system on a U. S. spacecraft. (Oleson, 
1999) | 

That particular Hall engine is a TAL (Thruster with Anode Layer) D-55 Hall- 
effect thruster. It was eked by the Russian Central Research Institute of Machine 
Building (TSNIIMASH), and it successfully operated, firing ten times from 23 to 24 
October. The engine provided 40 mN of thrust and raised the satellite’s orbit 0.35 NM. 
(Oleson, 1999) | 

As electrical propulsion ree NASA Glenn has clearly defined its objectives. 

For fiscal year 1999, four major product targets have been defined. Those targets 


include a next generation ion system for near- and further-term space science and a 
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pulsed plasma thruster (PPT) technology for precision imaging. In addition, NASA 
would like to develop both a high specific impulse, low-mass Hall system and a 
miniature electrothermal thruster for Earth science. (Oleson, 1999) 

NASA Glenn has also defined what it calls “grand challenges” in the development 
of electrical propulsion. Those challenges are delineated for each of the three types of 
electrical propulsion systems, and they are outlined as follows: (Oleson, 1999) 

Bicchothennal Thrusters: 

e integrated electronics / thruster package for microssececen applications 

Electrostatic Thrusters: — 

e lightweight, high total impulse ion systems for space science 

e very high total impulse ion system for far-term space science 

Electromagnetic Thrusters: 

e low-cost, simple, low-power Hall for miniature Earth-orbital spacecraft 

e high specific impulse (I,), direct drive Hall for space science and Earth orbital 

e increase thruster efficiency by factor of two to 16% 

e increase capacitor life to 40,000,000 pulses and beyond 

e lower electromagnetic emissions to below present standard specifications 

e predict contamination on spacecraft surfaces in design phase 

Nonetheless, the ability and timeliness of electric propulsion to meet these 
challenges is unpredictable. This is due inherently to the fact that the “state of the art” in 
electric propulsion is due and determined largely by the “state of the art” of the electrical 
components that comprise various aspects of the electrical propulsion system. Primarily, 


the concern is the required additional mass of electrical components. This fact will 
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become evident in the forthcoming discussion of the specific power or “alpha” (a) 
parameter. 
B. SUMMARY OF AVAILABLE ENGINE DATA 

One of the contributions of this thesis is to obtain and summarize the most current 
reported data regarding present-day electric propulsion engine performance 
characteristics. This summary, which includes data from various vendors as well as the 
above research sites, can be found on page 36 in Table 5: Summary of Current 
Technology aahine and the various sources of this information are listed at the bottom 
of the table. Figure 8, on page 35, provides an overall map of current engines and their 
respective values of specific impulse, I;. 

The electrical engine parameter that has certainly been the most difficult to define 
has been an accurate representation of the parameter a, specific power or “state of the 
art.” This parameter is so named because it defines the relation of the power conversion 
unit with respect to the system’s mass. In other words, as previously defined, 

a. = P. / mp, , where P, is the power supplied to the electrical propulsion system and mpp 
is the power-plant mass (Sutton, 1992, pg. 597). 

The apparent controversy surrounding the “state of the art” parameter is 
evidenced in the article “Optimization of Electric Propulsion Systems Considering 
Specific Power as Function of Specific Impulse” by M. Auweter-Kurtz et al. This article, 
found in the Nov.-Dec. 1988 issue of the AIAA Journal of Propulsion and Power, 
presents om argument that specific power, o, 1s a strong function of specific impulse, I;, 
or exhaust velocity, v. Specific impulse and exhaust velocity are related via the equation, 


v =I, go. 
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The article notes that several different organizations have conducted 
optimizations for various missions and that those optimizations assumed a constant 
specific power of the propulsion system. This assumption is arguably a potential source 
of error in that it implies the following about variations in power levels and specific 
impulse. First, it assumes that constant thrust and system efficiency could be evident 
with a changing specific impulse. Additionally, it implies that constant specific power 
could be maintained with changing power levels during optimization. (Auweter-Kurtz et. 
al., 1988, pg. 512) 

In undertaking this optimization analysis, this point was certainly considered. 
However, due to the fact that the trace of an engine’s efficiency is not continuous, 
average values were utilized over limited specific impulse ranges thereby assuming no 
strong dependence on specific impulse or exhaust velocity. 

Still, it is not elementary to assign « values to a particular engine. 

Items of concern in determining specific power include such things as the number 
of engines that will be necessary for a specific mission and whether or not those same 
engines can share various components of the power conditioning equipment. If power 
gonditenitg components can be shared, then « improves as mass is decreased. 
Likewise, a would be detrimentally affected if more mass were required such as 
additional solar cell panels to power required thrusters. 

Specific power also can be improved if the electrical propulsion system can be 
powered through the systems already integral to the payload. As such, the question arises 
as to where that mass should be billed. If it is attributed to the payload mass (m,)), then 


again the value for specific power is improved. 
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In summary of the above concerns, it should be evident that the a parameter is not 
only engine specific but also configuration (number of engines and shared components “i 
redundancy) and design (payload compatible power sources) specific. 

In order to demonstrate this aspect of the determination of the value for specific 
power, the following example is illustrated. 

Atlantic Research Corporation is one of the current suppliers of one of the 
commercially available Hall-type thrusters. In their brochure for the SPT-100, they 
depict the configuration for a “typical SPT-100 propulsion system for GEO satellites.” 

| The system depicted consists of the following major components: four SPT-100 
thrusters, four pairs of xenon flow controllers (XFC’s), two PPU-100 power processing 
units, the propellant management assembly (PMA), and two xenon storage tanks. The 
entire configuration allows for redundancy; therefore, we would assume that the designed 
mission would require the use of only two of the four thrusters. (ARC brochure, 1998) 

In determining specific power, it is necessary to account for the masses of all 
components that encompass the power-plant mass (Mpp). This does not include the actual 
mass of the propellant (m,). Thus, it is often referred to as the “dry mass” -- probably 
hailing from the early days of liquid rocket propulsion. 


The masses of all of the components are as follows: 


SPT-100 (4)..sccssessesseesssesessees 3.5kg/ SPT 
XEC (4 pairs)...scessssessseees0.64 kg pate 
PPU-100 (2)..ccscssessesssesseeeee. 6.2 kg / PPU 
PMA.....secssesseee Stee os 3.25 kg 


However, the mass of the xenon tanks is mission dependent. (Clauss, 1999) 
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The tank mass can be estimated in the manner demonstrated in Chapter IT, C. 
Again, assuming that a tank which holds zero propellant mass weighs 1 kg and that the — 
tank has a mass fraction of 0.1 to 0.12 (as a function of the mass of xenon), then for a 40 
kg of xenon, a 5 kg tank jen be required as 1 kg + (0.1)(40 kg) = 5 kg. (Clauss, 1999) 

In the brochure example for a GEO satellite, we pick a usable example of a 15- 
year GEO satellite for NSSK. Such a mission would require approximately 120 kg of 
xenon to be held in two tanks. Hence, 60 kg of xenon would be held in each tank. 
Accordingly, one tank would then weigh 7 kg as 1 kg + (0.1)(60 kg) = 7 kg. 

Nowhere in any current publications is there a discussion as to how to “charge” 
mass values in determining specific power or “alpha.” In this case, all “dead” mass is 
attributed to power-plant mass (mpp), and all mass that is proportional to power is 
assioned to payload mass (m,)). Hence, we would have to include 2 kg of tankage in our 
value for power-plant mass. 

With the value for tank mass, it is now possible to calculate the power-plant mass 
as follows: 

Mpp = 4 (3.5 kg)spr + 4 (0.64 kg)xFc + 2 (6.2 kg)ppy + 3.25 kg pma + 2.0 kg tanks. 
Thus, mpp = 34.2 kg. | 

Each of the SPT-100 thrusters requires 1.35 kW of power delivered at 300 V and 
4.5 A. With a known PPU efficiency of nppy = 0.93, we are able to calculate that a total 
of 5,808 W must be delivered to the thruster bank. Hence, we are now able to calculate 
the specific power for the configuration. (Clauss, 1999) 


a = P. / Mpp = 5,808 W / 34.2 kg = 169.8 W/kg. 
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It can not be overemphasized that this calculated value of specific power for the 
SPT-100 is only for the illustrated four-thruster configuration as designed for the mission 
of a 15-year GEO satellite. Otherwise, the masses change, and specific power must be 


re-calculated. 
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Figure 8: Specific Impulse (I,) Ranges for Current Technology Engines 
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Table 5: Summary of Current Technology Engines, 
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IV. PLANNING / OPTIMIZATION ALGORITHM — MATCHING THE 


MISSION, PAYLOAD, AND ENGINE 
OPTIMUM DESIGN PROCEDURE: 


1. Specify the Mission. 

Define the mission profile in terms of the velocity requirements (Au) and payload 
mass (mp). List additional constraints such as mission time (tm), power (Pe), cost, etc. 
2. Estimate the payload mass ratio, mp: / mo. 

Based upon past operational and flight-test experience, an initial estimate of 
suitable candidate electric propulsion thrusters may yield an approximate expected 
payload mass ratio. Note: This initial “ballpark” estimate is primarily based upon 
mission profile (velocity requirement) and suitable engines, looking at cost, mission time, 
and operational flight experience. _ 

‘In order to remain near the optimum for engine performance (as illustrated in 
earlier plots), a potential payload mass ratio range (from mp; = 0.35 to 0.55) is selected 
with the poeenondmne numbers for Au / v, and v/ Ve. 


IDEAL: without “tankage” penalty (i.e., particularly for Teflon PPT): 


Mpt / Mp = 0.35: Au /v, = 0.404; v/v, = 0.782. 
Mp1 / Mo = 0.45: Au/ Vv, = 0.327; v/v. = 0.820. 
Mp! / My = 0.55: Au/ Vv, = 0.257; v/v. = 0.864. 


“TANKAGE”: with 10 % tankage penalty: 


Mp1 / Mp = 0.35: Au/v,= 0.375; v/v. = 0.850. 
Mp1 / Mp = 0.45: Au/v, = 0.306; v/v, = 0.895. 
Mp / My = 0.55: Au / v, = 0.240; v/v. = 0.945. 
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If none of the listed payload ratios will suffice (or if further iteration shows that 
an intermediate value for payload ratio is sppropane): the required numbers for Au / v, 
and v / v, can be obtained from the plots in Figures 8 and 9. These graphs are for ideal 
systems and systems requiring “tankage” respectively 
3. Calculate the corresponding specific impulse, I,, based upon payload mass 
ratio... using the tabular data below and the equation, 

I, =v/ Qo. 


IDEAL: without “tankage” penalty: 


Mp; / My = 0.35: I, = 1.936 Au / Qo. 
Mp / My = 0.55: I, = 3.36 Au/ go, where gy = 9.81 m/s”. 


“TANKAGE”: with 10 % tankage penalty: 


Mp) / Mp = 0.35: I, = 2.27 Au/ Qo. 
Mp / Mo = 0.45: I, = 2.92 Au/ go. 
Mp / My = 0.55: I, = 3.94 Au/ go, where g, = 9.81 m/s”. 


For intermediate values of payload mass ratio, this same calculation can be 
accomplished using the ratios from Figures 8 and 9 along with the following equation, 

I, =(V/Ve)(Ve/ Au) (Au/ gp ). 

4, Designate a candidate electric propulsion engine. 

Referencing Figure 8 (on page 35), select the engine (or engines) whose available 
range of specific impulse, I,, most closely matches the calculated value necessary for the 
optimum profile. 

Engine selection also leeks in the characteristic engine performance values for 
specific power (a), thruster efficiency (n), a available thrust. These values can be 


found in the comprehensive summary of Table 5 on pages 36 and 37. 
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Again, if no engine proves an adequate match for the calculated values of I,, it is 
necessary to use the graphs depicted in Figures 9 and 10 and attempt intermediate values 
of mass payload ratio within the range. 

5. Calculate burn time, tp. 

Determine the thruster burn time for the profile from the following formula, 

tp = (( vc / Au ) Au)’ /(2 © Nt). 

At this point in the process, all of the above values are given: Au from mission 
profile, v, / Au from step 2, « and 7; from Table 5. 

6. Check that burn time, tp, is less than both mission time, t,,, as well as the design 
lifetime of the designated engine. 

If burn time eeesds acceptable mission time or design life of the thruster, select 
another electric propulsion candidate (step 4) or vary the payload mass ratio (step 2). 

7. Calculate the anticipated masses of the propellant (m,) and power-plant (mp,) via 
the following formulas for the respective payload mass ratios: 


IDEAL: without.“tankage” penalty (i.e., for Teflon PPT): 


Mp / Mp = 0.35: Mp = 1.152 mp1; Mpp = 0.612 m, 
Mp! / Mp = 0.45: Mp = 0.731 mp} ; Mpp = 0.672 my 
Mp / Mp = 0.55: Mp = 0.469 mp! 3 Mpp = 0.746 my 


“TANKAGE”: with 10 % tankage penalty (1.e., for all other engines): 


Mp! / Mp = 0.45: Mp = 0.679 mp1 ; Mpp = 0.801 mp 
Mpt / Mp = 0.55: Mp = 0.432 mp1 ; Mpp = 0.893 mp 


For intermediate payload ratios, the following equations must be used... 


(a) Mp1 / My = desired ratio. 


(b) Mp =(Vv/ve)* mp. 
(Cc) Mpo=M)+ Mppt+ Mp. 
(d) Mp) is given. 
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8. Determine the total mass of the assembled vehicle (m,) and verify the payload 
mass ratio. 

Mp = Mp + Mpp + MpI. mass ratio = Mp / Mo. 

9. Determine the required thrust based on the values for specific impulse, I,, and 
propellant mass, mp. Determine the required power to be delivered to the thruster. 

F =(mp)/tp ) I; go. Pe =Q Mpp. 

Note: This calculated value of P, should not greatly exceed the current “state of 
the art” for power generation (fuel cells, battery, or solar cells) of approximately 20 kW. 
If the power requirements cannot be met, return to step 4. 

10. From the thrust calculation and the engine thrust characteristics (step 4), 
calculate the number of engines required to provide that thrust for the assigned 
mission. 

11. Design check. 

Considering component redundancy, is this design acceptable? That is, with 
current technology, is this a reasonable abe of required engines? If not, return to 
design step 4 or attempt an intermediate payload ene (step 2). 

If yes, it may be advisable to iterate in order to account for redundancy thereby 
improving the a parameter. In other words, are there any components that the thrusters 
may be able to share, thereby reducing mp, and improving a. This will improve overall 
performance by reducing mass and possibly reducing the number of engines required for 


the mission. 


42 











0.8 ae 
viv, 
068 
° oe, 
= ee 
s Toe een. 
: ao 
0. 4 ° te ee Om 
ee ie ' , 
rea ~~ du 


anne OE Ee ilies ge Te ned iteeeen panes ene een momen 
[C3 0.36 0.36038 04 042 O44 0.48 048 05 O52 054055 058 OG 


Mpj / Mo 


Figure 9: Comprehensive Plot of Au/v, and v/v, for Ideal Systems. 


This plot serves as a baseline plot of the values of Au/ ve and v/v, for all 
anticipated payload mass ratios between the optimum selection points of 
~ Mp / my = 0.35 and 0.55. The points roughly bracket the optimum point at 
Mpi / Mo = 0.45. | 

This plot is usable only for ideal estimates of the calculations discussed in the 
design calculation, and it is usable for systems that do not require fuel tanks — i.e. no 
“tankage penalty.” Such electric thrusters or engines include the pulsed MPD thrusters: 


Teflon PPT, LES 8/9 PPT, NASA Primex EO-1, and the Primex PRS-101. 
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Figure 10: Comprehensive Plot of Au/v, and v/v. 
for Systems Requiring “Tankage” of 10%. 


This plot is acceptable for use with all electric propulsion mses that require 
fuel storage tanks. theicioe. it may be used for “tankage penalty” calculations for all of 
the engines listed except the pulsed MPD type. 

Because of the fact that these plots take fuel tank considerations into account, this 
chart (or at jens the determined values for the defined payload rations) should be 


consulted. 
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Figure 11: FLOWCHART -- Optimization Procedure for Electric Propulsion Engines 
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m, = 1.152 m,; 1.078m, 
Mp. = 0.612 m,; 0.723 m 


m, = 0.731 m,; 0.679 m, 2p: | 
Moe = 0.672 m,; 0.801 m, 


m =M+M,~+M, 
mass ratio = Ma /m, 


m=m+m,+m,y Mm, =m, +M, +My, 


mass ratio =m, /m, mass ratio=m, /m, 









F =(m,/t,) 1. ge 


F =(m,/t,)l g F =(m,/t,) | go 
R =a m,, 


P=am, | R=aKm, 





no m 
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| # engines = F/T 
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V. APPLICATION OF METHOD 


Note: In the following three examples, the defined parameters which had to be estimated 
by the author are preceded by the symbol “~”. 
A. “LIGHT” LEO TO GEO MISSION 
For a particular LEO to GEO satellite, the following parameters are defined: 
mission velocity requirement, Au = 4200 m/s 
payload mass, mp) = 250 kg 
mission time, t,, = unrestricted / reasonable. 
Using the defined optimization procedure, three potentially optimized engine runs 
are attempted at payload mass ratios of mp) / mM, = 0.35, 0.45, and 0.55. 
In this example, the second (“tankage” considered) case for the payload mass 
ratio of 0.45, proves to be the optimized profile. Note that approximately 10 to 12 
engines are required if the H2 arcjet is selected as the thruster for this profile. 
This example demonstrates that fuel tank mass considerations lead to the addition 
of more thrusters in order to sonnei for the added mass. The mass penalty is easily 


scalable. 
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. IDEAL: 


I, = 829 sec 

engine selected: 

NHs3 Arcjet 
I, = 650 sec 
a = 320 W/ke 
Tt = 0.32 
T=0.2N 


tp= 5.28 x 10° sec 
= 6.11 days 
ij< te * 


mp = 270 kg 

Mpp = 195.2 kg 

m, = 715 kg 
payload ratio checks 


F=3.26N; Pp =63 kW 
.17 engines required; 
power req. excessive 


with TANKAGE... 
I, = 972 sec 
engine selected: 
H> Arcjet 
I, = 1000 sec 
a = 333 
Me 0.4 
T=0.2-—0.25N 


tp= 4.71 x 10° sec 
= 5.45 days 


m, = 270 kg 

Mpp = 195.2 kg 

m, = 715 kg 
payload ratio checks 


F=5.47N; Pe=65 kW 
”.22 to 28 engines © 
required 


I, = 1075 sec 
engine selected: 
NASA / Primex EO-1 
| I, = 1150 sec 
a = 20 W/kg 
Nt = 0.098 
T = pulsed 


tp = 4.21 x 10’ sec 
= 487 days 
t:<ta 


“Mp = 182.8 kg 


Mpp = 122.8 kg 
My = 556 kg 
payload ratio checks. 


F = 0.0490 N; P,. = 3 kW 
. 1,919 engines 
required!!! 


I, = 1250 sec 
engine selected: 
H> Arcjet 
I, = 1000 sec 
a = 333 
Nt = 0.4 
T =0.2 —0.25N 


tp= 7.07 x 10° sec 
= 8.18 days 
tp < tm 


m, = 169.8 kg 
Mpp = 136.0 kg 
My, = 556 kg 


- payload ratio checks 


F=2.36N; Pe=46kW 
.10 to 12 engines 
required; power high 
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Mp / m= 0.55 


I, =.1439 sec 

engine selected: 

ARC SPT-100 
I, = 1600 sec 
a=~ 169.8W/kg 
At = 0.48 
T = 0.083 N 


tp = 1.638 x 10° sec 
= 19.0 days 
ty < tn 


m, = 117.3 kg 

Mpp = 87.5 kg 

My = 455 kg 
payload ratio checks 


F= 1.011 N; Pe=15 kW 
..13 engines required 


I, = 1687 sec 
engine selected: 
Hall (Xe) 
I, = 15-1600 sec 
a = 150 
t= 0.5 
T =0.04N 


tp = 2.04 x 10° sec 
= 23.6 days 
to<tn 


my = 108.0 kg 

Mpp = 96.4 kg 

m, = 454 kg 
payload ratio checks 


F=0.831N; P= 15 kW 


”. 21 engines required 











B. “HEAVY” LEO TO GEO (COMM SAT) MISSION 

Given a LEO to GEO mission for the orbit raising of a communication satellite, 
the optimization algorithm requires that we initially define the payload and required Au 
for the transit. | 

The payload is defined to be relatively massive, weighing 1000 kg, and the 
defined Au for the transfer is approximately 4640 m/s. For this mission, the transfer time 
cannot be unrealistically excessive due to the general mission urgency of communication 
satellite tasking. Using the planning algorithm, the results are illustrated as shown on the 
following page. 

This example provides some interesting results. First, the payload is extremely 
heavy for most electrical propulsion engines. As stated in the earlier chapters, most high 
thrust electric propulsion engines are stil in development. Yet, the Hall TAL D-55 
engine is capable of accomplishing this mission in 88% days with only 11 engines — 
excluding redundancy. 

In both of the cases at the payload _— ratio of 0.45, it is evident that the 
increased burn time to maintain the profile allows for fewer engines as less thrust is 
required over a longer time. None-the-less, the number of engines required for this 


profile as mp; / mp, = 0.45 is unacceptable. 
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IDEAL: 


I, = 916 sec 
engine selected: 
NASA/Primex EO-1 
s = 1150 sec 
a=~20 W/kg 
Nt = 0.098 
T = pulsed 


tp = 3.37 x 10’ sec 
= 389 days 
= 1.067 years 


SHOWSTOPPER... 
due tp >> tm 


with TANKAGE... 


I, = 1074 sec 
engine selected: 
H> Arcjet 
I, = 1000 sec 
a = 333 Wikg 
Nt = 0.4 
T =0.2-0.25N 


tp = 5.75 x 10° sec | 
= 6.65 days 
th 


m, = 1078 kg 
Mpp = 779 kg 
My = 2857 kg 
payload ratio checks 


F= 18.39 N; P= 259 kW 
..74 to 92 engines 
required!!! 


Mp! / My = 0.45 


I, = 1187 sec 

engine selected: 

Hall (Xe) 
I, = 15 -1600 sec 
a = 150 W/kg 


tp = 1.342 x 10° sec 
= 15.54 days 


Mp = 731 kg 

Mpp = 491 kg 

My = 2222 kg 
payload ratio checks 


F=8.02N; Pp=74kW 
..201 engines required! 


I, = 1381 sec 

engine selected: 

Hall (Xe) 
I, = 15 -1600 sec 
a = 150 W/kg 
Nit = 0.5 
T =0.04N 


tp= 1.533 x 10° sec 
= 17.74 days 


My = 679 kg 

Mpp = 544 kg 

My = 2223 kg 
payload ratio checks 


F=6.52N; Pp =82 kW 
.. 163 engines required! 
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Mpi / My = 0.55 


I, = 1589 sec 

engine selected: 

Hall TAL D-55 
I, = 1600 sec 
a=~ 50.9 W/kg 
| cee 0.48 
T = 0.082 N 


tp = 6.67 x 10° sec 
= 77.1 days 
(<1 


Mp = 469 kg 

Mpp = 350 kg 

m, = 1819 kg 
payload ratio checks 


F=1.104N; P.=18 kW 


. 14 engines required 


I, = 1864 sec 
engine selected: 
Hall TAL D-55 
I, = 1600 sec 
=~ 50.9 W/kg 
Nr = 0.48 
T = 0.082 N 


tp = 7.65 x 10° sec 
= 88.5 days 
its 


Mp = 432 kg 

Mpp = 386 kg 

My = 1818 kg 
payload ratio checks 


F = 0.886 N;. P= 20 kW 


_ ..11 engines required 








C. LEO TO MARS MISSION 

For a representative LEO to Mars mission, the optimization algorithm requires 
that we initially define the payload and required Au. 

The payload is defined to be a relatively small probe (mass equal to 250 kg), and 
the defined Au for the transfer is approximately 5,700 m/s. In this case, Au also defines a 
mission time (tm) equal to 256 days. Thus, in order for the calculations to be valid, the 
burn time (t,) must not exceed tm. 

The results of the calculations are demonstrated as shown. 

Note that for the third case of mp) / mM, = 0.55 either of the steady MPD thrusters 
would have been a better specific impulse fit. Those engines are scaleable in the range 
from I, = 2000 to 5000 sec. However, since they are presently still under development, 
they have been excluded from this analysis even though it is strongly anticipated that 


such engines will be the desired propulsion systems for deep-space missions. 
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IDEAL: 


I, = 1125 sec 

engine selected: 

Teflon PPT (pulsed) 
I, = 1000 sec 
a = 1.0 W/.kg 
Tit = 0.07 
T= pulsed 


tp= 1.422 x 10” sec 
= 16,457 days 
— = 45 years!!! 


SHOWSTOPPER.... 
due ty >> tm 


with TANKAGE... 


I, = 1319 sec 

engine selected: 

Hall (Xe) 
I, = 15-1600 sec 
a= 150 W/kg 
Nt = 0.5 
T =0.04N 


tp= 1.540 x 10° sec 
= 17.83 days 
Ath 


m, = 270 kg 

Mpp = 195.2 kg 

m, = 715 kg 
payload ratio checks 


F=2.58N; Pe=30kW 
.65 engines required; 
power required is high 


Mp1 / Mp = 0.45 


I, = 1458 sec 

engine selected: 

Hall SPT (Xe) 
I, = 1500 sec 
a= 150 W/kg 
Nt = 0.5 
T=0.04N 


tp = 2.03 x 10° sec 
= 23.4 days 
iota 


m, = 182.7 kg 

Mpp = 122.8 kg 

M, = 556 kg 
payload ratio checks 


F = 1.324N; Pp=19kW 
.34 engines required 


I; = 1697 sec 

engine selected: 

ARC SPT-100 
I, = 1600 sec 
a = 131.4 W/kg 
T = 0.083 N 


tp = 2.75 x 10° sec 
= 31.8 days 
tp <tmn 


My = 169.8 kg 

Mpp = 136.0 kg 

m, = 556 kg 
payload ratio checks 


F = 0.969 N; Pp = 18 kW 
.12 engines required 
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I, = 1952 sec 
engine selected: 
SPT-100 (Fakel) 
I, = 1600 sec 
a = 150 W/kg 
Nt = 0.48 
T = 0.083 N 


tp = 3.42 x 10° sec 
= 39.5 days 
ip ty 


m, = 117.3 kg 

Mpp = 87.5 kg 

My = 455 kg 
payload ratio checks 


F=0.538 N; Pe =14kW 


”.7 engines required 


I, = 2289 sec 

engine selected: 

XIPS 
I, = 28-3500 sec 
a = 100 W/kg 
t= 0.75 
T = 0.015-0.04 N 


tp = 3.76 x 10° sec 
= 43.5 days 
tp < tm 


m, = 108.0 kg 

Mpp = 96.4 kg 

m, = 454 kg 
payload ratio checks 


F = 0.789 N; Pp = 10 kW 
.20 to 53 engines 
required 











VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Routine implementation of electric propulsion for space thrusting. applications is 
relatively new even though the governing concepts are almost 40 years old. The very 
latest engines, presented as examples in this thesis, demonstrate that some aspects of 
electric propulsion are still in their early stages of development. However, in light of 
time-tested theory and the limitations of chemical propulsion systems, it 1s evident that, 
right now, electric propulsion is the only potential stepping-stone to deep-space probes or 
other craft that could conceivably travel the great distances between the planets and 
beyond the solar system. 

For any such mission profiles — interplanetary or interstellar in particular — the 
propulsion systems for these missions will have to be optimized. Such engines will have 
to be cost-effective, thrust efficient and therefore mass minimized. Someday there may 
even be considerations to make these systems potentially re-useable and refuelable. This 
thesis presents one technique with which mission profiles and engine selection can be 
matched and consequently optimized to theoretically achieve the greatest benefit of least 
overall mass with both maximum fuel and maximum payload in order to achieve the 
shortest mission time for a given payload. 

Thus, a primary focus of this thesis has been to provide a new way of looking at 
electrical propulsion systems design 1n addition to attempting to provide an up-to-date 
and comprehensive summary of available electric propulsion thruster data. The selected 
approach (i.e., the Langmuir-Irving formulation) is based upon first principles as well as 


an analysis of other methods that are currently in use. In essence, an objective has been 
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to produce a design alternative to the satellite design engineers. This methodology or 
“procedure” is intended to generate some consideration into performance near a pre- 
determined yet not inflexible “dual-optimum” point. It is apparent in the literature that 
aspects of the discussed “dual-optimum” point have been plotted in various texts and 
somewhat addressed but never fully explored. 

The most important soneept to note is the idea that these results are “theoretical” 
in that systems are currently designed according to algorithms that are based more solidly 
on operational experience and power considerations than on any approach that considers 
pure optimization. In essence, no design algorithm is really in use that seems to consider 
optimization. Granted, electric propulsion is still in its early stages and most applications 
involving satellite re-positioning and orbit-lifting do not place heavy lift nor aggressive 
profile constraints upon engines. Consequently, such considerations or additional 
restraints may not really be required. 

Due to the relatively limited number of electrical propulsion engines that are 
currently operational, the examples of this thesis show that it is sometimes difficult to 
find an identical match to an optimum profile’s specific impulse, I;. In fact, some of the 
most powerful and most dude electric propulsion engines — those most likely to be - 
desirable for high thrust and high payload mass missions — are still in some stages of 
development. 

Specific impulse, I,, is comparable to the required “gas mileage” for a particular 
trip through space. After the profile is determined through orbital mechanics, the | 
required incremental change in vehicular velocity (Au) is known. Then, the required tank 


of fuel can be established for an expected level of performance or fuel efficiency — “gas 
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mileage.” Hence, the objective is to match the profile to specific impulse (or engine). 
This is accomplished using the “dual-optimum” procedure. Nonetheless, the “dual- 
optimum” establishes not only I, but also suggests desired values of specific power (a) 
and thruster efficiency (n:) as well as the likely burn time (tp) for a desired mission 


profile, Au. 

As demonstrated in the three provided ee in the previous chapter, it is Bt 
obvious how a particular engine will perform on a specific mission profile. In essence, 
the design optimum is dependent as much on the payload mass (m,)) and mass payload 
ratio (Mp / M,) as upon the specific characteristic of the engine itself. 

The most important of these engine characteristics are specific power (a), thruster 
efficiency (7), and available thrust (T). Specific power, a, is the most critical to the 
designer because it drives the amount of available power to the thruster versus the 
required mass of the overall system that provides that power. im power conditioning 

systems improve and the “state-of-the-art” in the field advances, this parameter will 
improve and enhance electrical propulsion system er In this thesis, it is 
evident that specific power has been the most difficult parameter to establish for the 
majority of the listed engines. A worked example demonstrates that this is primarily due 
to the fact that the value given to « varies depending somewhat upon arbitrary 
considerations such as thruster configuration, redundancy, etc. In addition, it shows that 
mass reduction : equally as important as where that mass is “billed” within the 
calculations. 

Thruster efficiency, rn; , generates a feel for the engine’s efficiency in terms of 


nozzle losses, off-axis losses, as well as other thruster inefficiencies. Generally, these 
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values vary widely depending upon the type of thruster under consideration. For 
example, the LES 8/9 PPT thruster has an efficiency of only 0.0068 to 0.009 while the 
N>H, resistojet turns in numbers in the 0.8 to 0.9 range. The reasons for this are 
discussed in depth. 

Available thrust, T, although extremely small on most electrical propulsion 
systems, gives a measure of the number of engines required for a necessary profile. This 
value can be a major player as demonstrated in the examples because it drives largely the 
additional issues of component redundancy and additional mass — which again affects 
alpha, a. 

Another consideration which became apparent in the course of this work is the 
capability to diiow for the required fuel tanks within the design and thus within the 
optimum profile. Fuel considerations are certainly mission or profile dependent. 
Consequently, this is a substantial capability for the designer to predict optimum 
performance with variable tank mass as per mission requirements. 

These results, as presented, are easily scalable depending upon the particular 
profile and its fuel requirements. Exploration of the topic of allowing for “tankage” 
created estimation possibilities that allowed for even more realistic comparisons between 
different electric propulsion systems. As evident in the overall engine parameter 
summary in Table 5, some engines (of tie pulsed MPD variety) require no fuel storage 
tanks at all. | 

If the engine summary shows one thing, it certainly lets one know that there are a 
host of possibilities out there being developed by aggressive companies. These systems 


use a variety of physical principles to create thrust. Some at greater expense of the other 
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relevant parameters -- ©, 7, and T. Consequently, the examples provided for the design 
algorithm (on the three different mission profiles) demonstrate that all engines are not 
created equal for various mission profiles F certainly not if one endeavors to fly an 
optimum profile. | 

Of course, industry may argue that there are important additional considerations 
on any profile. For example, a communication satellite that takes six months to be in 
position is probably not cost-effective. That could be a mission failure. Nonetheless, 
success depends on what parameter is being optimized. In the case of this research, the 
objective was to optimize both the mass and time constraints. At the present time, in the 
development of the “state-of-the-art” in the evolution of electrical propulsion, this was 
not always possible. This fact is certainly demonstrated in the examples provided, which 
show that there are no obvious trends with the systems currently available. Electric 
propulsion is not totally there yet; however, scientists and engineers are pushing the 
frontier forward and expanding to new horizons and possibilities. 
B. RECOMMENDATIONS 

In any branch of research involving a rapidly growing technology (particularly 
one in which corporate vendors on the “cutting-edge” of that technology are especially 
unwilling to _ up data that they consider proprietary), it is extremely difficult to track 
the leading edge of developments. Consequently, acquiring the latest available data can 
be next to impossible. Hence, best-estimate scenarios can be flawed or only approximate. 

In the course of this research, the author frequently ran into barriers when 
attempting to answer questions regarding particular thruster systems. A more open forum 


could certainly benefit not only science but also the vendors themselves, particularly if 
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the data contained within this thesis proves useful to designers. On the other hand, the 
developed procedure might demonstrate that certain engines are less effective for | | 
necessary mission profiles. As a result, one design might prove less effective than 
another. 

Additional work could be extremely advantageous if a more in-depth study could 
be made into attempting to combine the various schools of thought in design criteria. 
Like missile design, there is no clear-cut or established method that is followed due to the 
fact that a design must consider the whole picture — 1.e., mass, time, and power | 
considerations —all at once. Considerations of power, time, and cost seem to be the 
biggest factors in design at the present time, and obviously contractors wanting to sell 
their products are not going to be quick to admit that a system could be better optimized 
with another company’s thruster. Certainly, perfection is infinitely expensive, but that is 
not what this process advocates. Systems could simply be better if any attempt were 
made to “ride the dual-optimum profile.” The author would certainly be curious if 


industry were to consider the presented ideas and improve upon them. 
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APPENDIX A. DEVELOPMENT - PAYLOAD MASS FRACTION EQUATION: 


By definition, specific power, & = Pe / Mpp. Thus, mpp = Pe / a. (Sutton, 1992, pg. 597) 


By definition, electrical power available (as previously shown), Pe= mp v’/(2 tp Nt). 
Hence, Mpp = (1/c) (mp v'/(2 tp N)). Solving for payload mass, mp, yields... 
Mp = Mpp (2 © ty Ne/ V’) OF Mpp/ Mp = V"/2 0 ty Mt. (Sutton, 1992, pg. 596) 


By definition of allocated spacecraft masses, mp = Mp + Mpi + Mpp. Define final mass, mg 
= Mp - Mp (i.e. initial mass less propellant mass). Thus, redefine initial mass as follows: 
My =M, + mz. (Sutton, 1992, pg. 596) 


Derived from the equation for maximum velocity at propellant burnout, 
Au =v In (m,/m)), know that e“”Y = m,/ mp=m, / (mo - mp)=1/(1-mp)/m,). This 
equation, when solved for mp / mo, gives Mp / Mp = 1 - ee“! (Sutton, 1992, pg. 123) 


Rewriting the equation for initial mass yields my = mp (1 + mpp / mp) + Mp). Now, when 
substituting the above equation for mpp / mp, the initial mass is rewritten as follows: 
Mp = Mp (1 +v/2a tp MH) + mp. 


Again, rewriting the equation... 1 = (mp / mp) (1 + v/2a tp Nt) + Mp1/ Mo. Substituting 
the equation for m, / m, and solving for mp; / mp, the equation transforms to 
Mp} / Mo = 1 -(1-e") (1 + v7 /2 & tp np). 


When this equation is multiplied through and solved for mpi/my , the payload mass 
fraction appears as illustrated: mp) / Mp = (1 - (v/ vy (eA - 1))/ eA”. 
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APPENDIX B. FIRST PRINCIPLES — CODES 


The following three pages contain the MATLAB codes that were utilized in the 


generation of Figures 2, 3, and 4 respectively. 
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% THESIS 

$ Plot2(2): mpl / mo vs. v/ve 

% "payload fraction vs. v / vc" 
% LT John Jay De Bellis, USN 

% 30 JUN 99 

$ variables: 

% x() =v / ve 

% a() = deltau / ve 

% y{) = mpl / mo 


xl=linspace (0,3, 300) 

x2=linspace (0, 3,300) 

x3=linspace (0,3, 300) 

al=0.2 

a2=0.3 

a3=0.5 

yl=(1-x1.*%2 .* (exp(al./x1)-1))./exp(al./x1) 
y2=(1-x2.*°2 .* (exp(a2./x2)-1))./exp(a2./x2) 
y3=(1-x3.*2 .* (exp(a3./x3)-1))./exp(a3./x3) 


figure(1), plot(xl, yl,*2;y2,x3;y3), grid 

axis([0 3 0.0 1.0]}) 

ylabel('mpl / mo'),xlabel('v / vc') 

gtext ('deltau/vce = 0.2'), gtext('deltau/vc = 0.3'), gtext('deltau/ve = 0.5") 
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% THESIS 
$ Plotl: deltau / ve vs. v / ve (or Is / ve) 
3 "incremental change of vehicle velocity vs. impulse" 
¢ LT John Jay Dé Bellis, USN 
$ 25 AUG 99 7 
Variables: 
x{) =v / ve 
a() = mpl / mo 
y() = deltau / ve 


A WP al ol 


xl=linspace(0,3,300) 
x2=Linspace (0,3,300) 
x3=linspace (0,3,300) 


al=0.35 
a2=0.45 
a3=0.55 


Vilexl:*log( (1-4 x1. °2).:7 Gal. -F- xis*2)) 
V2=x2 3" LOG (1. R282) f Lad. R222) 
V5=x3.* log (¢(1.-4+--x3..°2).2/ (a3. + -43.72)) 


Ligure (1). 7 plottxl, yi, xZ7 V2; 3, v3); orid 
axis({0.0 2.0-0.0°0.5]) 


ylabel('deltau / vc'),xlabel('v / ve') 
gtext ('mpl/mo = 0.35'), gtext('mpl/mo = 0.45'), gtext('mpl/mo 
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0. 55)") 


% THESIS 

$ Plot3: optimization curves -- joint optimum : 

% "product of incremental change of vehicle velocity and payload ratio vs. 
% specific impulse" (mpl/mo * deltau vs. Is) 

%$ LT John Jay De Bellis, USN 

$ 25 AUG 99 

$ Variables: 

% x() =v/ve 

gz a() = mpl / mo 

% y() = deltau / ve 


x1l=linspace (0,3,300) 
x2=linspace (0, 3,300) 
x3=linspace (0,3, 300) 


al=0.35 
a2=0.45 
a3=0.55 


yl=xl.*log({1 + x1.*2)./(al + x1l.%2)) 
V2Z=x2."* 10g CL +: 22672) 07 (a2 + 222°2):) 
Vo=x3e* Log ( (1. -* x3.°2) 27 las + 32°27} 


$zl=yl.*al 
$z2=y2.*az 
$z3=y3.%*a3 


figure (1); plotixzl, yl, x2; y2,%37 3) +. 9rad 

axe ClOcl 0.15025 1.5 )) 

ylabel('deltau/vc * mpl/mo'), xlabel('v/vc') 

gtext ('mpl/mo = 0.35'), gtext('mpl/mo = 0.45'), gtext('mpl/mo = 0.55") 
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APPENDIX C. DUAL OPTIMUM CODE 


This appendix lists the necessary Maple program code that is utilized in the 


production of Figure 4, which illustrates the “dual optimum.” 
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THESIS: Plot #5: OPT. CURVES: 3-D plot: optimization of v/vc and mpl/mo curves 


LT John Jay De Bellis, USN 
30 JUN 1999 


defined variables: 
deltau = defined incremental change of vehicle velocity (m/s) -- defined by mission 


vc = characteristic velocity (m/s) 
v = effective exhaust velocity (m/s) 
mratio = payload ratio = mpl/mo 


X = v/vc 
y = mpl/mo 


> restart; with(plots); with(plottools); 
plot3d(x*y*1n((14x%2) / (y+x42)),x=0.5..1.3,y=0.3..0.6,axes=BOXED, ti 
tle="optimization: product of v/vc and mpl/mo"); 

[ animate, animate3d, animatecurve, changecoords, complexplot, complexplot3d, conformal, 
contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, display, display3d, 
fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, listcontplot, 
listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logpiot, matrixplot, odeplot, pareto, 
pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, 
polyhedraplot, replot, rootlocus, semilogplot, setoptions, setoptions3d, spacecurve, 
sparsematrixplot, sphereplot, surfdata, textplot, textplot3d, tubeplot | 

[ arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, 
ellipticArc, hemisphere, hexahedron, homothety, hyperbola, icosahedron, line, octahedron, 
pieslice, point, polygon, project, rectangle, reflect, rotate, scale, semitorus, sphere, stellate, 


tetrahedron, torus, transform, translate, vrml] 
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APPENDIX D. TANKAGE PENALTY CODES 


The following pages list the MAPLE program codes that were used in the creation 


of Figures 6 and 7 respectively. 
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THESIS 
Plot1T: deltau / ve vs. v / ve (or Is / vc) with 10% TANKAGE PENALTY 


"incremental change of vehicle velocity vs. impulse” 


OP AO oO - 


LT John Jay De Bellis, USN 


A JO 


25 AUG 99 
$ Variables: 
% x() =v /ve 
% a() = mpl / mo 
% y() = deltau / ve 


xl=linspace (0,3, 300) 
x2=linspace (0,3, 300) 
x3=linspace (0,3, 300) 


al=0.35 
a2=0.45 
a3=0.55 


yl=xl.*log((1.1 + STAD) iS (al + x1.2 % O45) 
V2=x2.* log ( (lak + e26°2) 67 (a2 x2.°2 + 0.1) 
Vo=xo Log (led XS. °2) 07 (asd: PS 2: 021)) 


figure(]):,- plot (x1, v1) 2; V2; S57 V3) 5 grid. 

axi-$( [0.0 32:0 0.0 :0.5)) 

ylabel('deltau / vc'),xlabel('v / vc') 

gtext('mpl/mo = 0.35'), gtext('mpl/mo = 0.45'), gtext('mpl/mo = 0.55') 
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% THESIS 
$ Plot1iT: deltau / ve vs. v / ve (or Is / ve) with/without 10% TANKAGE PENALTY 


% "incremental change of vehicle velocity vs. impulse” 


%* LT John Jay De Bellis, USN 
% 25 AUG 99 


$ Variables: 

% x() =v /ve 

3 a() = mpl / mo 

% y() = deltau / ve 


xl=linspace (0,3, 300) 
x2=linspace (0,3,300) | 
x3=linspace (0,3,300) 


a=0.45 


Vix le LOG (lisa RAS 2) ey (as SORA, SOE 7) 
y2=x2.*log((1 + x2.%2)./(a + x2.%2)) 


figure(1), plot(xl,yl,x2,y2), grid 

axis{(0.0°2.0° 0.0 0.51) 

ylabel('deltau / vc'),xlabel('v / ve') 

gtext('mpl/mo = 0.45"'), gtext('mpl/mo = 0.45 with 10% tankage penalty') 
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